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t:;;j- . Abstract 

We describe some field theoretic methods for studying quantum spin systems in 
^ \ one dimension. These include the nonlinear cr-model approach which is particularly 

H ' useful for large values of the spin, the idea of Luttinger liquids and bosonization 

^ , which are more useful for small values of spin such as spin- 1/2, and the technique of 

j^ \ low-energy effective Hamiltonians which can be useful if the system under considera- 

tion is perturbatively close to an exactly solvable model. We apply these techniques 
, - to similar spin models, such as spin chains with dimerization and frustration, and 

^ , spin ladders in the presence of a magnetic field. This comparative study illustrates 

O ' the relative strengths of the different methods. 

> 

(N 

oo . 

^ ■ 1 Introduction 

o. 

One-dimensional and quasi-one-dimensional quantum spin systems have been studied ex- 
tensively in recent years for several reasons. Many such systems have been realized experi- 

^ \ mentally, and a variety of theoretical techniques, both analytical and numerical, are avail- 

able to study the relevant models. Due to large quantum fluctuations in low dimensions, 
'^ ■ such systems often have unusual properties such as a gap between a singlet ground state 

and the excited nonsinglet states; this leads to a magnetic susceptibility which vanishes 

O ■ exponentially at low temperatures. Perhaps the most famous example of this is the Hal- 

dane gap which was predicted theoretically in integer spin Heisenberg antiferromagnetic 
'k>( \ chains |1[, and then observed experimentally in a spin-1 system Ni{C2lI^N2)2N02{ClOi) 

^ ' 0. Other examples include the spin ladder systems in which a small number of one- 

dimensional spin- 1/2 chains interact amongst each other ||^. It has been observed that if 
the number of chains is even, i.e., if each rung of the ladder (which is the unit cell for the 
system) contains an even number of spin-1/2 sites, then the system effectively behaves 
like an integer spin chain with a gap in the low-energy spectrum. Some two-chain ladders 
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which show a gap are {VO)2P207 i, SrCu203 g and Cu2{C^Hi2N2)2Ch 1- Con- 
versely, a three-chain ladder which effectively behaves like a half-odd-integer spin chain 
and does not exhibit a gap is Sr2CuzO^ 0. A related observation is that some quasi- 
one-dimensional systems such as CuGeO^ spontaneously dimerize below a spin-Peierls 
transition temperature [0; then the unit cell contains two spin-1/2 sites and the system 
is gapped. Another interesting class of systems are the alternating spin chains such as 
bimetallic molecular magnets. An example is NiCu{pbaOH){H20)s-2H20 in which spin- 
I's {NP~^) and spin-l/2's {Cu^~^) alternate. The ground state of these systems have a 
nonzero total spin 5*0. It turns out that there is a gap to states with spin greater than So, 
but no gap to states with spin less than 5*0. 

The results for gaps quoted above are all in the absence of an external magnetic field. 
The situation becomes even more interesting in the presence of a magnetic field |^ . Then 
it is possible for an integer spin chain to be gapless and a half-odd-integer spin chain to 
show a gap above the ground state for appropriate values of the field ^, |TU|, |TT], |T2|, |T^ . 
This has been demonstrated in several models using a variety of methods such as exact 
diagonalization of small systems, bosonization and conformal field theory |T^, [T^, and 



perturbation theory [|T6]. In particular, it has been shown that the magnetization of 
some systems can exhibit plateaus at certain nonzero values for some finite ranges of the 
magnetic field. 

The plan of this paper is as follows. In Sec. 2, we discuss the low-energy properties of 
the dimerized and frustrated antiferromagnetic spin chain. In Sees. 3 and 4, we present 
some field theoretic methods which can be used for studying spin chains and ladders with 



or without an external magnetic field ^7j. These methods rely on the idea that the low- 
energy and long- wavelength modes of a system (i.e., wavelengths much longer than the 
lattice spacing a if the system is defined on a lattice at the microscopic level) can often 
be described by a continuum field theory. In Sec. 3, we discuss the nonlinear a-model 
approach, while in Sec. 4, we discuss the concepts of Tomonaga-Luttinger liquids and 
bosonization. In Sec. 5, we discuss the low-energy effective Hamiltonian approach and 
show how it can be combined with bosonization to gain an understanding of the magnetic 
properties of one-dimensional spin systems. 



2 Spin Chain with Dimerization and Frustration 

Experimental studies of some of the quasi-one-dimensional spin systems have shown that 
besides the nearest neighbor antiferromagnetic exchange, there also exists a second neigh- 
bor exchange J2 of the same sign and comparable magnitude. Such a second neighbor 
interaction has the effect of frustrating the spin alignment favored by the nearest neighbor 
interaction. Therefore, a realistic study of one-dimensional systems requires a model with 
both frustration (J2) and dimerization (governed by a parameter 6). 

The Hamiltonian for the frustrated and dimerized antiferromagnetic spin chain can be 



written as 



H 



^ [ 1 - (-1)^ 5]S,- S^,+i + J2Y. S^■ S,+2 , 



where the hmits of the summation depend on the boundary condition (open or periodic). 
(We have set the average nearest neighbor interaction Ji to be equal to 1 for convenience). 
The interactions are schematically shown in Fig. 1. The region of interest is defined by 
J2 > and < (5 < 1. 




Figure 1: Schematic picture of the spin chain described by Eq. ([1|). 

The ground state properties of the Hamiltonian (|l|) have been studied at some rep- 
resentative points in the J2 ~ S plane using the density-matrix renormalization group 



(DMRG) method [18|. The phase diagrams obtained for spin-1/2 and spin-1 chains are 
shown in Figs. 2 and 3 [|r^. We use the word 'phase' only for convenience to distin- 
guish between regions with different modulations of the two-spin correlation function as 
discussed later. Our model actually has no phase transition even at zero temperature. 

For the spin-1/2 chain I^U], ^, the system is found to be gapless on the line A which 
runs from J2 = to J2C = 0.241 for 6 = 0; see Fig. 2. The model is gapped everywhere 
else in the J2 — S plane. There is a disorder line B given hj 2J2 + 6 = 1 on which the exact 
ground state of the model is given by a product of singlets formed by the nearest-neighbor 
spins which are joined by the stronger bonds {1 + 6); this is called the Shastry-Sutherland 
line [^, and it ends at the Majumdar-Ghosh point (J2 = 0.5,6 = 0). The correlation 
length ^ goes through a minimum on B. Finally, the peak in the structure factor S{q) is 
at Qmax = vr to the left of B (called region I), decreases from vr to 7r/2 as one goes from B 
up to the line C (region II), and is at qmax = 7r/2 to the right of C (region III). 

In the spin-1 case (Fig. 3), the phase diagram is more complex. There is a solid line 
marked A which runs from (0, 0.25) to about (0.22 ± 0.02, 0.20 ± 0.02) shown by a cross. 
To within numerical accuracy, the gap is zero on this line and the correlation length ^ is 
as large as the system size N. The rest of the 'phase' diagram is gapped. However the 
gapped portion can be divided into different regions characterized by other interesting 
features. On the dotted lines marked B, the gap is finite. Although ^ goes through a 




Figure 2: Ground state phase diagram of the spin-1/2 chain in the J2 — S plane. 



maximum when we cross B in going from region II to region I or from region III to region 
IV, its value is much smaller than N. There is a dashed line C extending from (0.65, 0.05) 
to about (0.73, 0) on which the gap appears to be zero (to numerical accuracy), and ^ is 
very large but not as large as A^. In regions II and III, the ground state for an open chain 
has a four-fold degeneracy (consisting of S" = and S = 1), whereas it is nondegenerate 
in regions I and IV with S = 0. The regions II and III, where the ground state of an 
open chain is four-fold degenerate, can be identified with the Haldane phase. The regions 
I and IV correspond to the non-Haldane singlet phase. Regions I and IV are separated 
by the disorder line D given by 2 J2 + 6 = 1, while regions II and III are separated by line 
E. The lines B, D and E seem to meet in a small region V where the ground state of the 
model is numerically very difficult to find. 

As can be seen from Fig. 1, setting 6 = 1 results in a two-chain ladder where the 
interchain coupling is 2 and the intrachain coupling is J2. We can hold J2 fixed and vary 
the interchain coupling J. Numerical studies show that for spin-1/2, the system is gapped 
for any nonzero value of J, although the gap vanishes linearly as J ^ 0; this can be shown 
using bosonization. On the other hand, the spin-1 chain has a finite value of the gap for 



any value of J |T^ 




Figure 3: Ground state phase diagram of the spin-1 chain. 



3 Nonlinear cr-model 



The nonhnear a-niodel (NLSM) analysis of spin chains with the inclusion of J2 and 6 
proceeds as follows p3||. We first do a classical analysis in the S* — > 00 limit to find the 
ground state configuration of the spins. Let us make the ansatz that the ground state is 
a coplanar configuration of the spins with the energy per spin being equal to 



eo 



S' 



(1 + 5) cos 6*1 



[l-S) cos^s + J2 cos(^i + 62) 



(2) 



where 6*1 is the angle between the spins S2J and S^2i+i and 6*2 is the angle between the 
spins S2i and 821-1- Minimization of the classical energy with respect to the 9i yields the 
following three phases. 

(i) Neel: This phase has 61 = 62 = vr; hence all the spins point along the same line and 
they go as TiTi along the chain. This phase is stable for 1 — 6^ > 4J2- 

(ii) Spiral: Here, the angles 61 and 62 are given by 



cos 6*1 



1 

TT5 



1-6^ 
4J2 



6 
1 + 52 



4J2 



and cos Oo 



1-5 



4J2 



5 



1-52 



4Jo 



(3) 



where 7r/2 < ^1 < tt and < ^2 < ^i- Thus the spins he on a plane. This phase is stable 
for 1 - 52 < 4J2 < (1 - 5^)1 5. 

(iii) Colinear: This phase (which needs both dimerization and frustration) is defined to 
have 9i = 71 and 62 = 0; hence all the spins point along the same line and they go as ttii 
along the chain. It is stable for (1 — 6'^)/S < 4J2. 

These phases along with their boundaries are depicted in Fig. 4. Thus even in the classical 
limit S ^ 00, the system has a rich ground state 'phase diagram'. 
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Figure 4: Classical ground state phase diagram of the spin chain with frustration and 
dimerization. 

We can go to the next order in 1/S, and study the spin wave spectrum about the 
ground state in each of the phases. The main results are as follows. In the Neel phase, we 
find two zero modes, i.e., modes for which the energy Uk vanishes linearly at certain values 
of the momentum k, with the slope duk/dk at those points being called the velocity. The 
two modes are found to have the same velocity in this phase. In the spiral phase, we 
have three zero modes, two with the same velocity describing out-of-plane fluctuations, 
and one with a higher velocity describing in-plane fluctuations. In the colinear phase, we 
get two zero modes with equal velocities just as in the Neel phase. The three phases also 
differ in the behavior of the spin-spin correlation function S{q) = J2n{So ■ Sn) exp(— ign) 
in the classical limit. S{q) is peaked at g = (^1 + ^2)/2, i.e., at g = vr in the Neel phase, at 
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7r/2 < g < 7r in the spiral phase and at g = 7r/2 in the cohnear phase. Even for S = 1/2 
and 1, DMRG studies have seen this feature of S{q) in the Neel and spiral phases W^ . 
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Figure 5: Classical configuration of the spins in the Neel phase. 

We now derive a NLSM field theory which can describe the low-energy and long- 
wavelength excitations. In the Neel phase, this is given by a 0(3) NLSM with a topological 
term [|l], |15|. The field variable is a unit vector which is defined as follows. The classical 
ground state in the Neel phase has a unit cell, labeled by an integer n, with two sites 
labeled as In and 2n respectively; see Fig. 5. We define linear combinations of the two 
spins as 

r S-in — S; 



L 



In — J2n 

2S 

Sin + S2n 

2a 



(4) 



Here a is the lattice spacing; hence, the size of each unit cell is 2a. Note that 

0, 



'n ' ^n 



s 52 ' 

so that (pn becomes an unit vector in the large S limit. These fields satisfy the commutation 
relations 



[I 



ma ) Ynb 



2a 



Omn / , f 



abc (Pnc i 



(6) 



where m, n are unit cell labels, a, 6, c denote the components x, y, z, and eabc is the com- 
pletely antisymmetric tensor with exyz = 1- This means that we can write In = <Pn^ n„, 

— # — * 

where the vector 11 is canonically conjugate to </>, i.e., 

[ 0ma , n„b ] = — 5mn ^ab ■ (7) 

We now go to the continuum limit by introducing a spatial coordinate x which is equal 
to 2na at the location of the rfi^ unit cell. Summations get replaced by integrals, i.e., 
Sn '~^ J dx/{2a). The commutation relation (|^) then takes the form 

[ M^) , My) ] = ^S{x-y) 6at . (8) 

We note that and 0' are orthogonal to because is an unit vector. We will see below 
that both / and 11 are given by first-order space-time derivatives of 0. In the low-energy 
and long- wavelength limit, the dominant terms in the Hamiltonian will be those which 
have second-order derivatives of 0, and therefore first-order derivatives of /. To find this 
Hamiltonian, we rewrite (0) in terms of and /, and Taylor expand these fields to the 
necessary order, i.e., 

0„+i = (pix) + 2a0'(x) + 2a20" + ... , 
r„+i = fix) + 2af{x) + ... , (9) 

where x = 2na. We then use the constraints in (^ and do some integration by parts 
(throwing away boundary terms at a; = ±cxd) to obtain the continuum Hamiltonian 

H = [ dx[^ (r+ ^0')' + 7^ 0" ] , (10) 

where 



/ 


dx 

c 
9 


eg' 
2 


2aS^/l 


$r + 


c 

2g' 




-4J2 

2 


-6^ 


and 


27r5(l 


4J2- 
-6). 


F' 



(11) 

By expanding ([lO|) to second order in small fluctuations around, say, = (0, 0, 1), we find 
an energy-momentum dispersion relation of the 'massless relativistic' form uj = c\k\; thus 
c is the spin wave velocity. Similarly, by expanding (|l^) to fourth and higher orders in 
small fluctuations, we find that g' is the coupling constant governing the strength of the 
interactions between the spin waves. 

One can show that the Hamiltonian (|10D follows from the Lagrangian density 
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(Incidentally, one can derive the canonically conjugate momentum 11 and then the angular 
momentum / from (p!2|), 



n = 


1 ^ 9 ^ ^, 

cg^ An 




r = 


X n = — - 6x6- 
cg^ 


9 
A-K 



(13) 

thereby verifying that / and 11 only contain first-order derivatives of as stated above). 
From (p!2D, we see that 9 is the coefficient of a topological term, because the integral of 
this term is an integer which defines the winding number of a field configuration 0(x, t)). 
For 9 = IT mod 27r and g less than a critical value gc, it is known that the system is gapless 
and is described by a conformal field theory with an SU{2) symmetry |T5|, |2^. For any 



other value of 9, the system is gapped, and the gap is of order AE ~ exp{—2n/g'^). For 
J2 = S = 0, one therefore expects that integer spin chains should have a gap of the order 
exp(— vrS") (note that this goes to zero rapidly as S ^ 00, so that there is no difference 
between integer and half-odd-integer spin chains in the classical limit), while half-integer 
spin chains should be gapless. For the two-spin equal-time correlation function, this 
means that < Sq ■ Sn > should decay as a power-law (— l)"/|ra| as \n\ -^ 00 for half-odd- 
integer spin chains, and exponentially as (— l)"exp(— ra/^) for integer spin chains, where 
the correlation length C, ~ c/AE. All this is known to be true even for small values of 
S like 1/2 (analytically) and 1 (numerically) although the field theory is only derived for 
large S. In the presence of dimerization, one expects a gapless system at certain special 
values of 6. For S = 1, the special value is predicted to be 6c = 0.5. We see that the 
existence of a gapless point is correctly predicted by the NLSM. However, according to 
the DMRG results, 6c is at 0.25 for J2 = |^ and it decreases with J2 as shown in Fig. 
3; this differs from the NLSM results in (|Tl|) according to which 9 should be independent 
of J2. These deviations from field theory are probably due to higher order corrections in 
1/5 which have not been studied analytically so far. 

In the spiral phase of the J2 — 6 model, it is necessary to use a different NLSM which 
is known for 6 = |^6|, ^. The field variable is now an 5*0(3) matrix R. The Lagrangian 



density is 

C = ^tr(E^EPo) - ^tr(E'^E'Pi), (14) 



where c = S{l + y)^/T^^/y, g"^ = 2y'(l + y)/{l - y)/S with l/y = AJ2 , and Pq and Pi 
are diagonal matrices with diagonal elements (1, 1, 2y{l — y)/{2y^ — 2y + 1)) and (1, 1, 0) 
respectively. Note that there is no topological term; indeed, no such term is possible since 
n2(50(3)) = unlike U2{S^) = Z for the 0(3) NLSM in the Neel phase. Hence there 
is no apparent difference between integer and half-integer spin chains in the spiral phase. 
A one-loop renormalization group [^ and large N analysis [^ indicate that the system 



should have a gap for all values of J2 and 5, and that there is no reason for a particularly 
small gap at any special value of J2 ■ The 'gapless' point found numerically at J2 = 0.73 
for spin-1 is therefore a surprise. 
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Finally, in the colinear phase of the J2 — 5 model, the NLSM is known for 6 = 1, i.e. 



for the spin ladder. The Lagrangian is the same as in (|T3), but with c = AaS J J2{J2 + 1), 
g^ = Jl + I/J2/S and ^ = 0. There is no topological term for any value of 5*, and the 



model is therefore gapped. 

The field theories for general 6 in both the spiral and colinear phases are still not 
known. Although the results are qualitatively expected to be similar to the 6 = case in 
the spiral phase and the 6=1 case in the colinear phase, quantitative features such as 
the dependence of the gap on the coupling strengths require the explicit form of the field 
theory. 

The NLSMs derived above can be expected to be accurate only for large values of the 
spin S. It is interesting to note that the 'phase' boundary between Neel and spiral for 
spin-1 is closer to the the classical {S — >■ 00) boundary 4J2 = 1 — 6"^ than for spin-1/2. 
For instance, the cross-over from Neel to spiral occurs, for 5 = 0, at J2 = 0.5 for spin-1/2, 
at 0.39 for spin-1, and at 0.25 classically. 

To summarize, we have studied a two-parameter 'phase' diagram for the ground state 
of isotropic antiferromagnetic spin-1/2 and spin-1 chains using the NLSM approach, and 
have compared the results with those obtained numerically. We find that the spin-1 dia- 
gram is considerably more complex than the corresponding spin-1/2 chain with surprising 
features like a 'gapless' point inside the spiral 'phase'; this point could be close to a critical 
point discussed earlier in the literature p^ ^. It would be interesting to establish this 
more definitively. 

Our results show that frustrated spin chains with small values of S exhibit some 
features not anticipated from large S field theories like the NLSMs. The NLSMs also 
leave many questions unanswered. For instance, the 0(3) NLSM which is applicable in 
the Neel phase does not tell us the exponent of the gap which opens up as one moves 
away from 9 = n (for g < gc) or as we go across g = gc (for 9 = n). To address these 
questions, we have to use the more powerful technique of bosonization. 

The NLSM approach can also be used to study spin chains in the presence of a magnetic 
field. Consider adding a Zeeman term to the Hamiltonian in (|I|), i.e., 

H = E [ 1 - i-^y ^]S^■ S,+i + J2Y. S^■S,+2 - J2 B-S,, (15) 

i i i 

where B denotes the magnetic field. In the region 1 — 5^ > 4 J2, the classical ground state 
of this Hamiltonian is given by a coplanar configuration in which the spins 821 and 5*21+1 
lie at angles ± a respectively with respect to the magnetic field, so that the angle between 
the spins 821 and 821+1 is 2a. Minimization of the energy fixes the angle a to be 

« = cos-^(^) • (16) 

(We are assuming that |-B| < 45, otherwise all the spins will align with the magnetic field 

10 



and a will be zero). We now define 

lb sm a 

In = ^1^^±^. (17) 

2a 

Note that the definition of (p is slightly different from the one in (^) in order to ensure that 
is an unit vector. However, / is orthogonal to and has the same commutation relations 
with as before. We can now go to the continuum limit and derive the Hamiltonian 



f (' ^ IJ"^' ^ 2? '" 



H = I dx[-S^(l + i^^'f + ^_4,'^-B-l 



where 



/ 



c = 2aS'sina yl — 4J2 — <^^ 
2 



Ssniayjl — 4J2 — 5"^ 
and e = 27rSsma (1-6) . (19) 

We can show that this follows from the Lagrangian density 

C = -1^ [0x<^+B-(B-0)0]2 - ^0'2 ^ A0.0'x<^, 

' ' +2Bx$-$+B'^ -{B-$f] ^ 0'2 ^ 0. 0' X . (20) 



We see from this that 

L[$x$ + B - (5-0)0] - / 



r= 0xn = — ^ [0x0 + 5 - (5-0)0] - — 0' . (21) 



Since c(7 = 4a and 5 ■ = in the classical ground state, we see that / is equal to 
B/{4a) plus small fluctuations; this agrees with its deflnition in (|1^ and the classical 
conflguration of the spins. One can now analyse the fleld theory governed by (|20|) using 
the renormalization group and other methods. We refer the reader to |^ for further 
details. 



4 Bosonization 

A very useful method for studying spin systems in one dimension is the technique of 
bosonization. Before describing this method, let us briefly present some background 
information. Further details can be found in Refs. |1^, |15|, Q plj] . 
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In one dimension, a great variety of interacting quantum systems (both fermionic 
and bosonic) is described by the Tomonaga-Luttinger hquid (TLL) theory. Typically, 
a TLL describes quantum systems which are translation invariant and gapless, i.e., the 
excitation energy above the ground state is zero in the limit of the system size L -^ 
oo. A TLL differs in three significant ways from the well-known Fermi liquid theory 
which describes many fermionic systems in two and three dimensions. First, all the 
low-energy excitations in a TLL have the character of sound modes which are bosonic 
and have a linear dispersion relation between the energy and the momentum (with the 
constant of proportionality being the sound velocity v). Even if the underlying theory is 
fermionic, the low-energy excitations are given by particle-hole pairs which are bosonic. 
The properties of a TLL are governed by two important parameters, namely, an interaction 
parameter K (noninteracting systems have K = 1) and the velocity v. Secondly, the one- 
particle momentum distribution function n{k) for fermions, which is obtained by Fourier 
transforming the fermion Green's function 



G{x,t) = <0|TV^(a;,t)V^f(o,0)|0> 



(22) 



and computing the residue of its pole in the complex u plane as a function of k, has no 
discontinuity at the Fermi surface k = kp for a TLL. Instead, it has a cusp there of the 
form 

n{k) = nikp) + const. sign{k - kp) \k - kpl^^'^^'^^^ . (23) 

On the other hand, in a Fermi liquid, n{k) has a finite discontinuity at the Fermi surface; 
see Fig. 6. Finally, correlation functions in a TLL typically decay at large distances as 
power-laws which depend on K, unlike the correlation functions of a Fermi liquid where 
the power-laws are universal. 




Figure 6: One-particle momentum distribution function for (a) an interacting Fermi liquid, 
and (b) a Tomonaga-Luttinger liquid. 

Let us be more specific about the nature of the low-energy excitations in a one- 
dimensional system of interacting fermions. Assume that we have a system of length 
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L with periodic boundary conditions; the translation invariance and the finite length 
make the one-particle momenta discrete. Suppose that the system has Nq particles with 
a ground state energy Eq{Nq) and a ground state momentum Pq = 0. We will be inter- 
ested in the thermodynamic limit A^o^ L ^ oo keeping the particle density po = Nq/L 
fixed. If we could switch off the interactions, the fermions would have two Fermi points, 
at A; = ±kF respectively, with all states with momenta lying between the two points being 
occupied. (See Fig. 7 for a typical picture of the momentum states of a lattice model 
without interactions). Even in the presence of interactions, it turns out that the low-lying 
excitations consist of two pieces P^ , 

(i) a set of bosonic excitations each of which can have either positive momentum q or 
negative momentum —q with an energy e^ = vq, where < q « kp and v is the sound 
velocity, and 

(ii) a certain number of particles Nji and N^ added to the right and left Fermi points 
respectively, where Nf{,NL « Nq. Note that Nfi and Nl can be positive, negative or 
zero. 

The quasiparticle excitations in (i) have an infinite number of degrees of freedom 
(in the thermodynamic limit), and they determine properties such as specific heat and 
susceptibility to various perturbations. The particle excitations in (ii) only have two 
degrees of freedom and therefore play no role in the thermodynamic properties. The 
Hamiltonian and momentum operators for a one- dimensional system (which may have 
interactions) have the general form 

H = EoiNo) + j:vq[ h^^^hn,, + h^^L,, ] 

+ MiVi? + iVL) + ^(^R + ^L? + ^(iV^-iVL)', 
P = E ? [ ~^i,Mi - bUh^ ] + [kp + jiNn + Nl) ] (Nn - Nl) , (24) 

where q is the momentum of the low-energy bosonic excitations created and annihilated 
by 6j and bg, K is a positive dimensionless number, and fi is the chemical potential of 
the system. We will see later that v and K are the two important parameters which 
determine all the low-energy properties of a system. Their values generally depend on 
both the strength of the interactions and the density. If the fermions are noninteracting, 
we have 

V = vf and K = 1 . (25) 

Note that one can numerically find the values of v and K by varying Nji and A^^^ and 
studying the 1/L dependence of energy and momentum of finite size systems. 

The technique of bosonization (combined with conformal field theory) is very useful 
for analytically studying a TLL |]T^ ^, ^, ^. This technique consists of mapping 
bosonic operators into fermionic ones, and then using whichever set of operators is easier 
to compute with. 



13 



To begin, let us consider a fermion with both right- and left-moving components. 
We introduce a chirality label u, such that u = R and L refer to right- and left-moving 
particles respectively. Sometimes we will use the numerical values z/ = 1 and —1 for R 
and L; this will be clear from the context. Then the second quantized Fermi fields are 
given by 



where Uk = 0, ±1, ±2, 



^Ax) = 


-1 oo 

V^ fc=_oo 


k = 


27r 


and 




{Cu,k, c^' 


,fc'} = , 


\Cu,k^ Cyi 


hi] = 5vv' Skk' ■ 



(26) 



(27) 



Next we define bosonic operators 



V-q fc=-oo 
-| oo 

Kq = -?=- E 4,fc-gC^,fc • (28) 

V'-q fc=-oo 

Note that 6^^ and 6^^ ^ create excitations with momenta q and —q respectively, where the 
label q is always taken to be positive. We can show that 

^v,q-, K',q'] = , and [K^g, bl,^^,] = S^^i 5,^/ . (29) 

The vacuum state of the system is defined to be the state |0 > which is annihilated by 
the operators c^^k for /c > and cli^ for /c < 0, and therefore by bu,q for all q. 

Let us define the chiral bosonic fields 

Mx) = -^ E ^ [ K, e^^''--""/^ - big e-'^--"''/^ ] _ v^ iv. , (30) 

■^V^ q>0 V^q ^ 

where the length parameter a is a cut-off which is required to ensure that the contribution 
from high-momentum modes do not produce divergences when computing correlation 
functions. The fields in (|30|) satisfy 



[0^(x),0i,/(x)] = - — (5,,,,/ sign (x - x) (31) 

in the limit a — > 0. It is useful to define two fields dual to each other 

0(x) = (I)r{x) + 0l(x) , 

0{x) = - (t)n{x) + 0l(x) . (32) 
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Then [(^{x) , (p^x')] = [9{x),9{x')] = 0, while 

[0(x),^(x')] = ^sign(a;-x'). (33) 

Now it can be shown that the ferniionic and bosonic operators discussed above are 
related to each other as 

^R = -^ VR e-''^'^'^ , 

i^L = -^VLC^'^f^, (34) 

in the sense that they produce the same state when they act on the vacuum state |0 >, 
and they have the same correlation functions. The unitary operators rjfi and rji are called 
Klein factors, and they are essential to ensure that the fermionic fields given in Eq. ( |5^ ) 
anticommute at two different spatial points x and y. 

The densities of the right- and left-moving fermions are given by p^ = ipji^pR ^-nd 
Pl = '4'l'^l- The total fermionic density and current are given by 

1 90 

P - Po = Pr + PL = - ^ TT ^ 

'TT OX 



J = vf{pr-Pl) = -T^Yx^ ^^^^ 



where po is the background density (fluctuations around this density are described by the 
fields ip or (p), and the velocity vp will be introduced below. 

Let us now introduce a Hamiltonian. We assume a linear dispersion relation 6,^ ^ = vpk 
for the fermions. The noninteracting Hamiltonian then takes the form 

oo 

Ho = VF J2 k[c^R^kCR,k + clkCL,k] + -^{Nl + Nl) 

k = —OD 

= -VF ("^ dx [iji{x)id^i,R{x) - ij{{x)id^Mx)] + ^{Nl + Nl) (36) 
JO L 



in the fermionic language, and 



^t . , J,t ,. N , ^^^, 



Ho = VFY.<i( ^rMs + bicpL,, ) + -^{Nl + Nl) 

q>0 ^ 

= VF r dx [ {d,<PRf + {d,<PLf ] 

= ^ r dx[ (5.0)2 + {djf ] (37) 

2 JO 

in the bosonic language. 
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We now study the effects of four-fermi interactions. Let us consider an interaction of 

tlie form 

1 fL 
V = - dx [2g2 Pr{x)pl{x) + g^ ( p|(x) + pi(x) ) ] . (38) 

Physically, we may expect an interaction such as gp^ jl-, so that g2 = 94, = 9- However, it 
is instructive to allow (72 to differ from g^ to see what happens. Also, we will not assume 
anything about the signs of (72 and (74. In the fermionic language, the interaction takes 
the form 



1 00 

^ = ^ Y. [ '^92C^RM+k3(^RMci,k2+k:,CLM 

ki,k2,k3=—oo 



+ 9i\'^R,kx+k3^RM^R,k2-ki^RM + ^LM+ka^LM^LM-ki^LMh 



(39) 

From this expression we see that 92 corresponds to a two-particle scattering involving both 
chiralities; in this model, we can call it either forward scattering or backward scattering 
since there is no way to distinguish between the two processes in the absence of some 
other quantum number such as spin. The g^ term corresponds to a scattering between 
two fermions with the same chirality, and therefore describes a forward scattering process. 

The quartic interaction in Eq. ( PUJ ) seems very difficult to analyze. However we 
will now see that it is easily solvable in the bosonic language; indeed this is one of the 
main motivations behind bosonization. The bosonic expression for the total Hamiltonian 
H = Hq + V is found to be 



92 



H = y: 9bF(0«,,+&L^Lj+|^(&kA.+^^A,.)+|^(^k.^i?.a+&iA,,)] 

q>0 ^^ ^^ 



+ ^{Nl + Nl)+ 'j-NnN, + ^{Nl + Nl). 



(40) 



The 9^ term only renormalizes the velocity. The g2 term can then be rediagonalized by a 
Bogoliubov transformation. We first define two parameters 



K 



, 9i 92 . , . 9i , 92 . 

271 2ti ' ^ 277 2ti ' 

, 9a 92 . , I , 9a , 92 

"^+2^-2^)/^"^+2^ + 2^ 



1/2 



1/2 



(41) 



Note that K < 1 if (72 is positive (repulsive interaction), and > 1 if (72 is negative (attractive 
interaction). [If (72 is so large that vp + 5'4/(27r) — g2/{2T:) < 0, then our analysis breaks 
down. The system does not remain a Luttinger liquid in that case, and is likely to go into a 
different phase such as a state with charge density order]. The Bogoliubov transformation 
then takes the form 



m,q 



jR,q 






L,q 
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7 bL,q + 7 4,, 

Kg = 



2 



1 - K 

where 7 = — , (42) 

' 1 + K ^ ^ 

for each value of the momentum q. The Hamihonian is then given by the quadratic 
expression 

(j>0 

77?? 1 -S /S A A 

+ ^[;^(iVH + iV,)2 + i^(iV«-iV^f ]. (43) 

Equivalent ly, 

H = \ 1"^ dx[ vKIi^ + |. (9,0)2 ] . (44) 

z JO A 

The old and new fields are related as 



(l + K)(t)R - (1- 


'K)cPl 


2VK 

(1 + iT) 0i - (1 - 


-K)4>n 



2VK 

K 6 and 9 = -^ 6 . (45) 



Note the important fact that the vacuum changes as a result of the interaction; the new 
vacuum |0) is the state annihilated by the operators hy^q. Since the various correlation 
functions must be calculated in this new vacuum, they will depend on the interaction 
through the parameters v and K. In particular, we will see below that the power-laws of 
the correlation functions are governed by K. 

Given the various Hamiltonians, it is easy to guess the forms of the corresponding La- 
grangians. For the noninteracting theory [g^ = fy'4 = 0), the Lagrangian density describes 
a massless Dirac fermion, 

£ = iil)R{dt + VFd^)i)R + iipL{dt-VFdr,)ipL (46) 

in the fermionic language, and a massless real scalar field. 



in the bosonic language. For the interacting theory in Eq. (]4^), we find from Eq. (145|) 
that 
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Although the dispersion relation is generally not linear for all the modes of a realistic 
system, it often happens that the low-energy and long-wavelength modes (and therefore 
the low-temperature properties) can be described by a TLL. For a fermionic system in 
one dimension, these modes are usually the ones lying close to the two Fermi points with 
momenta ^kp respectively; see Fig. 7. Although the fermionic field ip generally has 
components with all possible momenta, one can define right- and left-moving fields ipji 
and ipL which vary slowly on the scale length a, 



i,{x,t) = i,R{x,t) e^'^^^ + i^L{x,t)e 



-ikpx 



(49) 



Quantities such as the density generally contain terms which vary slowly as well as terms 
varying rapidly on the scale of a, 



P- Po 



i/j'ip 



^Ur + i^^L + e-'2/c,x ^y^ ^ ^2k,. ^y^ 



IT dx 



2Txa 



[r7jjr7ie*(2v/¥<^-2fe^x) ^ ^^^^g-i(20F<^-2fc^x) ]_ ^^q^ 




Figure 7: Picture of the ground state of a one-dimensional system of noninteracting 
fermions on a lattice. Filled circles denote occupied states lying below the Fermi energy 
Ef = 0. 



One can now compute various correlation functions in the bosonic language. Consider 
an operator of the exponential form 



O. 



^i2y^{m<f>+ne) 



(51) 
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(Such an operator can arise from a product of several ip^s and ip^^s if we ignore the Klein 
factors; then Eq. ( P^ implies that m ± n must take integer values). We then find the 
following result for the two-point correlation function at space-time separations which are 
much larger than the microscopic lattice spacing a, 

{0\TO^,n{x,t)Oly{0,0) |0) 

nn ^^^_^ _ i(^^[g-^[l^Ymy/K-n/yfK?[yt + X- iasign(t))(™^/^+"/^/^)' 

(52) 

Note that the correlation function decays as a power-law, and the power depends on 
the interaction parameter K. In the language of the renormalization group, the scaling 
dimension of Om,n is given by 

do = m^K + — . (53) 

K 

We can now discuss a spin chain from the point of view of bosonization. To be specific, 
let us consider a spin-1/2 chain described by the anisotropic Hamiltonian 

^ = E [ ^ (S^Sr,, + S-St,) + JAS^SU, - hSn, (54) 

where the interactions are only between nearest neighbor spins, and J > 0. S^ = S^ +iSf 
and S~ = Sf — iSf are the spin raising and lowering operators, and h denotes a magnetic 
field. Note that the model has a U{1) invariance, namely, rotations about the 5*^ axis. 
When A = 1 and /i = 0, the U{1) invariance is enhanced to an SU{2) invariance, because 
at this point the model can be written simply as H = J Y,i Si ■ S^+i. 

Eq. ([5^ ) is the well-studied XXZ spin-1/2 chain in a longitudinal magnetic field. 
It can be exactly solved using the Bethe ansatz, and a lot of information can then be 
obtained using conformal field theory [0, Q. The following results are relevant for us. 
The model is gapless for a certain range of values of A and h/J. For instance, this is 
true if — 1 < A < 1 and ^ = 0; then the two-spin equal-time correlations have oscillatory 
pieces which decay asymptotically as 



{Sq d„ ) 

WoSn) 



1)" 



_l)n 



\n 



1 1 

where r] = - + — sin^^ (A) . (55) 

2 vr 

For A > 1 and h = 0, the system is gapped; there are two degenerate ground states which 



have a period of two sites consistent with the condition (|68|). Thus the invariance of the 
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Hamiltonian under a translation by one site is spontaneously broken in the ground states. 

This is particularly obvious for A — i> oo where the two ground states are H 1 — ■ ■ ■ and 

— I \- ■ ■ ■. The two-spin correlations decay exponentially for A > 1 and h = 0. Finally, 

the system is gapped for h/J > 1 + A with all sites having S^ = 1/2 in the ground state, 
and for h/J < — 1 — A with all sites having Sz = —1/2. 

However, it is not easy to compute explicit correlation functions using the Bethe 
ansatz. We will therefore use bosonization to study the model in (|5^) . 

We first use the Jordan- Wigner transformation to map the spin model to a model of 
spinless fermions. We map an | spin or a | spin at any site to the presence or absence of 
a fermion at that site. We introduce a fermion annihilation operator ipi at each site, and 
write the spin at the site as 

S^ = VlV. - 1/2 = n, - 1/2 

Sr = (-i)^^,e^-E,". , (56) 

where the sum runs from one boundary of the chain up to the {i — l)*'^ site (we assume 
an open boundary condition here for convenience), nj = or 1 is the fermion occupation 
number at site i, and the expression for S^ is obtained by taking the hermitian conjugate 
of S'j", The string factor in the definition of S~ is added in order to ensure the correct 
statistics for different sites; the fermion operators at different sites anticommute, whereas 
the spin operators commute. 

We now find that 

H = - J2[i i4^^+l + h.c.) + J A (n,-l/2)(n,+i-l/2) - h{n,- 1/2)] . (57) 
i ^ 

We see that the spin-flip operators S^ lead to hopping terms in the fermion Hamiltonian, 
whereas the S^S'^j^^ term leads to an interaction between fermions on adjacent sites. 

Let us first consider the noninteracting case given by A = 0. By Fourier transforming 
the fermions, ipk = J2j '4^je~'^^^°' /\/N , where a is the lattice spacing and the momentum k 
lies in the first Brillouin zone — vr/a < k < vr/a, we find that the Hamiltonian is given by 

H = J2 ^k 4^k , (58) 

k 

where 

ujk = —J cos{ka) — h . (59) 

The noninteracting ground state is the one in which all the single-particle states with 
Uk < are occupied, and all the states with Uk > are empty. If we set the magnetic 
field h = 0, the magnetization per site rn = J2i Sf/N will be zero in the ground state; 
equivalently, in the fermionic language, the ground state is precisely half-filled. Thus, for 
m = 0, the Fermi points {uk = 0) lie at ka = ±7r/2 = kpa. Let us now add the magnetic 
field term. In the fermionic language, this is equivalent to adding a chemical potential 
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term (which couples to rti or Sf). In that case, the ground state no longer has m = 
and the fermion model is no longer half-filled. The Fermi points are then given by -^kp-, 
where 

kpa = TT (m + -) . (60) 

It turns out that this relation between kp (which governs the oscillations in the correlation 
functions as discussed below) and the magnetization m continues to hold even if we turn 
on the interaction J A, although the simple picture of the ground state (with states filled 
below some energy and empty above some energy) no longer holds in that case. 

In the linearized approximation, the modes near the two Fermi points have the veloci- 
ties dujk/dk = ±f , where v is some function of J, A and h. Next, we introduce the slowly 
varying fermionic fields ijjr and ipi as indicated above; these are functions of a coordinate 
X which must be an integer multiple of a. Now we bosonize these fields. The spin fields 
can be written in terms of either the fermionic or the bosonic fields. For instance, S^ is 
given by the fermion density as in Eq. (BB) which then has a bosonized form given in Eq. 



(0). Similarly, 

S+{x,t) = (-!)"> [V^Jj(x,t)e-''^^"/" + V^l(x,t)e''=^"/"] X 

where (—1)^/" = ±1 since x/a is an integer. This can now be written entirely in the 
bosonic language; the term in the exponential is given by 

1 fx 1 

dx'ip''{x',t)%jj{x',t) = — / dx'dx'(j) = — —j=^ [0ii:(a^, i^) + 0L(a^, i^)] , (62) 



where we have ignored the contribution from the lower limit at x' = — oo. 

We can now use this bosonic expressions to compute the various two-spin correlation 
functions G''\x,t) =< 0|T^'^(x, t)^^(0, 0)|0 >. We find that 



C{x,t) = m^ + ci 
G+-(x,t)+G~+(x,t) = C3 



+ 



cos(2A;i;'x) 

+ C2- 



{x + vty {x-vtyi "(x^-vH^)^ 

t_l\x/a 



+C4 



(X2 - Z;2t2)l/4K 



+ 



^,2 _ ^2^2-\(i7^-v^)^ 4X — ft)2 {x + Vt) 



(63) 



where ci, ...,C4 are some constants. The Luttinger parameters K and v are functions of 
A and h/J (or m). [The exact dependence can be found from the web site given in Ref. 
||1CI|| ; this contains a calculator which finds the values oi R = l/y/AnK and h/J if one 
inputs the values of M = 2m and A]. For h = 0, K is given by the analytical expression 

^ = 1 + -sin-i(A) . (64) 

K TT 
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Note that at the SU{2) invariant point A = 1 and h = 0, we have K = 1/2, and the two 
correlations G^^ and G^^ have the same forms. 

In addition to providing a convenient way of computing correlation functions, bosoniza- 
tion also allows us to study the effects of small perturbations which may take the system 
away from a TLL. For instance, a physically important perturbation is a dimerizing term 

V = SY: i-iy [ ^ {SrS._,, + Si St,,) + J^StSl, ] , (65) 

where 8 is the strength of the perturbation. Upon bosonizing, we find that the scaling 
dimension of this term is K. Hence it is relevant if A' < 2; in that case, it produces an 
energy gap in the system which scales with 5 as 

AE ~ 5^/(2-^) . (66) 

For the isotropic case A = 1, we have i^ = 1/2 and the gap scales as 1S.E ~ S^l"^ . [This 
is the exponent of the gap which appears as we vary 5 to move away from the gapless 
line (0 < J2 ^ <^2c; (^ = 0) for spin-1/2 in Fig. 2 or the line A for spin-1 in Fig. 3]. 
This phenomenon occurs in spin-Peierls systems such as CuGeO^\ below a transition 
temperature Tgp., they go into a dimerized phase which has a gap. 

Another interesting perturbation occurs when the frustration parameter J2 crosses the 
critical value J2C = 0.241 for 5 = in the spin-1/2 chain; see Fig. 2. This turns out to 
be a marginal perturbation, and it produces a gap which has an essential singularity of 



the form AE ~ exp[— const. /(J2 — J2c)] |33]. Because of this form, it is very hard to 



numerically measure the gap if J2 is close to J2c- 

Finally, when two isotropic spin-1/2 chains (with the spin variables in the two chains 
being denoted by S"^^^ and S^'^) are coupled together with a weak interchain coupling 

^ = ^'E5i^^-^^ (67) 

n 

we find that the perturbation S^^^ ■ S^'^^ has the scaling dimension 1. Hence this perturba- 
tion is relevant, and it produces an energy gap which scales as AE ~ J' . This has been 



confirmed by numerical calculations [jT9|. 



5 Low-energy Effective Hamiltonian approach 

As mentioned in Sec. 1, a quantum spin system can sometimes exhibit magnetization 
plateaus. For a Hamiltonian which is invariant under translation by one unit cell, the 
value of the magnetization per unit cell is quantized to be a rational number at each 
plateau. The necessary (but not sufficient) condition for the magnetization quantization 
is given as follows [^. Let us assume that the magnetic field points along the z axis, the 
total Hamiltonian H is invariant under spin rotations about that axis, and the maximum 
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possible spin in each unit cell of the Hamiltonian is given by S. Consider a state ip such 
that the expectation value of 5*2 per unit cell is equal to nis in that state, and ip has a 
period n, i.e., it is invariant only under translation by a number of unit cells equal to n or 
a multiple of n. (It is clear that if n > 2, then there must be n such states with the same 
energy, since H is invariant under a translation by one unit cell). Then the quantization 
condition says that a magnetic plateau is possible, i.e., there is a range of values of the 
external field for which ip is the ground state and is separated by a finite gap from states 
with slightly higher or lower values of total Sz, only if 

n { S — nig ) = a.ia integer. (68) 

Note that the saturated state in which all spins point along the magnetic field trivially 
satisfies (^) since it has rris = S (or —S) and n = 1. 

In this section, we study the magnetization as a function of the applied field for a two- 
and three-chain ladder using a perturbatively derived low-energy effective Hamiltonian 
(LEH) [jT^, ^ . In both cases, the first-order LEH will turn out to be the model described 



in Eq. (Q). As we pointed out earlier, a lot is known about this model [|10| ^. In 
particular, we will see that the exponent t] for the correlation power laws can be read off 
from the expression for the first-order LEH. 

We consider a three-chain spin-1/2 ladder governed by the Hamiltonian 

3 ^ ^ 3 

H = J' "^ Yl ^"-^^ ' Sa+l,n + J Y Y Sa,n " Sa,n+l " ^ X! X! ^l,n i (69) 

a n a=l n a=l n 

where a denotes the chain index, n denotes the rung index, h denotes the magnetic field, 
and J, J' > 0; see Fig. 8. We may choose h > since the region h < can be deduced 
from it by reflection about h = 0. It is convenient to scale out the parameter J, and 
quote all results in terms of the two dimensionless quantities J' /J and h/J. We will only 
consider an open boundary condition in the rung direction, namely, the summation over 
a in the first term of ( pUD runs over 1, 2. 

We now discuss the LEH approach for studying the properties of spin ladders. There 
are two possible limits which may be considered. One could examine J' /J -^ which cor- 
responds to weakly interacting chains, and then directly use techniques from bosonization 
and conformal field theory; this has been done in detail by others ||10|, |lT|. We therefore 
consider the strong-coupling limit J/ J' -^ which corresponds to almost decoupled rungs. 
In that limit, the LEH has been derived to first order in J/ J' for a three-chain ladder with 



periodic boundary condition along the rungs fSg, ^ , and for a two-chain ladder [jl^, |34 

We derive the LEH as follows. We first set the intrachain coupling J = and consider 
which of the states of a single rung are degenerate in energy in the presence of a magnetic 
field. In general, there will be several values of the field, denoted by h^, for which two 
or more of the rung states will be degenerate ground states. We will consider each such 
value of ho in turn. The degenerate rung states will constitute our low-energy states. 
If the degeneracy in each rung is d, the total number of low-energy states in a system 
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Figure 8: Schematic picture of the three-chain ladder described in Eq. 
1, 2 and 3 denote the three chains. 



. The labels 



with L rungs is given by L°'. (In general, the number d depends both on the system and 
on the field h^. It is two for both the models we will study here). Next, we decompose 
the Hamiltonian of the total system a.s H = Hq ^- V , where Hq contains only the rung 
interaction J' and the field /zq, and V contains the small interactions J and the residual 
magnetic field h — ho which are both assumed to be much smaller than J' . Let us denote 
the degenerate and low-energy states of the system as Pi and the high-energy states as 
Qa- The low-energy states all have energy Eq, while the high-energy states have energies 
Ea according to the exactly solvable Hamiltonian Hq. Then the first-order LEH is given, 
up to an additive constant, by degenerate perturbation theory. 



H. 



(1) 



eff 



J2 \Pi) {Pi\^\Pj) (Pj 



(70) 



V 



The calculation of the various matrix elements in Eqs. (^) can be simplified by using 
the symmetries of the perturbation V, e.g., translations and rotations about the z axis. 

To derive the LEH for the three-chain ladder, we decompose the Hamiltonian in (|69|) 
as H = Ho + V , where 



Ho 
V 



a=l,2 n a=l n 

3 ^ ^ 3 

J Z2 Yl '^a.n ■ Sa,n+l " (^ " ^o) zJ zJ ^a,n 



(71) 



a=l n 



a=l n 



We determine the field ho by considering the rung Hamiltonian Ho and identifying the 
values of the magnetic field ho where two or more of the rung states become degenerate. 
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The eight states in each rung are described by specifying the S^ components (+ and 
— denoting +1/2 and —1/2 respectively) of the sites belonging to chains 1, 2 and 3. For 
instance, the four states with total S = 3/2 are denoted by |1), ..., |4), where |1) = | + ++) 
and the other three states can be obtained by acting on it successively with the operator 
S^ = J2a ^a ■ These four states have the energy J'/2 in the absence of a magnetic field. 

There is one doublet of states |5) and |6) with S = 1/2, where |5) = [ 2 | H h) — | — 

++) — I + H — ) ]/\/6 and |6) ~ 'S'~|5). These have energy —J'. Finally, there is another 
doublet of states |7) = [ | + H — ) — | — 1-+) ]/v^ and |8) ~ S~\7) which have zero energy. 
It is now evident that the state |1) with S^ = 3/2 and the state |5) with 5*^ = 1/2 become 
degenerate at a magnetic field ho = 3J'/2, while states |5) and |6) are trivially degenerate 
for the field ho = 0. We now examine these two cases separately. 

For ho = 3J'/2, the low-energy states in each rung are given by |1) and |5), while the 
other six are high-energy states. We thus have an effective spin-1/2 object on each rung n. 
We introduce three spin-1/2 operators {S^, S^, S^) for each rung such that S^ = S^±iS^ 
and S^ have the following actions: 

S~ |1)„ = |5)„ , S- |5)„ = , 
and ^^|1)„ = ^|1)„, S:|5)„ = - ^|5). . (72) 

Note that the state which has a |1) on every rung, i.e., |111 ■ ■ ■), is just the state with 
rung magnetization rris = 3/2 corresponding to the saturation plateau. The state with a 
1 5) on every rung corresponds to the rris = 1/2 magnetization plateau. The LEH we are 
trying to derive will therefore describe the transition between these two plateaus. 

We now turn on the perturbation V in ([7T| ) with the assumption that J and h — ho 



are both much smaller than J'. We can write V = Y.n Ki,n+i, where 

3 ^ ^ 1 3 

Vn,n+1 = J 2^ Sa,n-Sa,n+1 — 7^{h — h^) 2_^ [ ^l^n + ^^l^n+l] ■ (73) 

a=l ^ a=l 

The action of Vn^n+i on the four low-energy states involving rungs n and n + 1 can be 
obtained after a long but straightforward calculation. We then use Eq. ([70|) and find that 
the LEH to first order in J/ J' is given, up to a constant, by 



2 



n n+1 



- ( /^ - ^ - ^ ) E ^^ (74) 

n 

where we have substituted ho = 3 J'/2. Thus the LEH up to this order is simply the XX Z 
model with anisotropy A = 1/2 in a magnetic field h — 3J'/2 — J/2; see Eq. (|5^) . 

We now use (|7^) to compute the values of the fields hi and /i2 where the states with 
all rungs equal to |1) and all rungs equal to |5) respectively become the ground states. 
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We can then identify hi with the lower critical field h^- for the plateau at rris = 3/2, and 
/i2 with the upper critical field hc+ for the plateau at rris = 1/2. 

To compute the field hi, we compare the energy Eq of the state with all rungs equal 
to |1) with the minimum energy Emin{k) of a spin- wave state in which one rung is equal 
to 1 5) and all the other rungs are equal to |1). A spin wave with momentum k is given by 

|fc) = 4^ E e^'" I5n) , (75) 

where |5n) denotes a state where only the rung n is equal to |5). The spin-wave dispersion, 
i.e., u{k) = E{k) — Eq, is found from (|7^ to be 

1 3 7' 7 

uik) = J{ cask - -^) + (h - — - -) . iJQ) 

This is minimum at A; = vr and it turns negative there for h < hi, where 

S 7' 

hi = — + 2J. (77) 

This is therefore the transition point between the ferromagnetic state |111---) and a 
spin-wave band lying immediately below it in energy. 

Similarly, we compute the field /i2 by comparing the energy Eq of the state with 
all rungs equal to |5) with the minimum energy Emin{k) of a spin wave in which a |5) 
at one rung is replaced by a |1). For a spin wave with momentum k, the dispersion 
uj{k) = E{k) — Eq is found to be 

1 72 9 ^ 3 7' 7 

uj(k) = J ( cosA; --) + —(--— C0S2A; ) - {h - — - -) . (78) 
^ ' ^ 2 ^ J' ^ 9 18 ' ^ 2 2 ' ^ ' 

This is minimum at A; = vr and it turns positive there for h > h2, where 

"^ T 

h, = — - J . (79) 

This marks the transition between the state |555 ■ ■ ■) and the spin- wave band. Equation 
([79|) agrees to this order with the higher-order series given in the literature |jlO| . 



From the first-order terms in (|7^, we can deduce the asymptotic form of the two- 
spin correlations. From (^5]), we see that the exponent r] = 2/3 for A = 1/2. Although 
this is the exponent for the H — correlation of the effective spin-1/2 defined on each 
rung, we would expect the same exponent to appear in all the correlations (5'^/5'^„) 
studied by DMRG in the previous section, regardless of how we choose the chain indices 
a,b = 1, 2, 3. We find that the analytically predicted exponent of 2/3 agrees quite well 



with the numerically obtained exponents which lie in the range 0.61 to 0.70 [^ . 

We now consider the LEH at the other magnetic field /iq = where the rung states 
1 5) and |6) are degenerate. We take these as the low-energy states and introduce new 
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effective spin-1/2 operators for each rung with actions similar to Eqs. (|72D, except that 
we replace |1) and |5) in those equations by |5) and |6). We compute the action of the 
perturbation V on the low-energy states, and deduce the LEH to be 

Heff = JY. Sn-Sn+1 - hJ2 S:,. (80) 

n n 

This Hamiltonian describes the transition between the magnetization plateaus at rus = 
1/2 and rris = —1/2; since these plateaus are reflections of each other about zero magnetic 
field, it is sufficient to study one of them. By a calculation similar to the one used to 
derive ([77|), the field hi can be found from the dispersion of a spin wave in which one 
rung is equal to |6) and all the other rungs are equal to |5). The dispersion is 

u{k) = h + J { cosk - I) . (81) 

This gives 

hi = 2J . (82) 

This is the lower critical field he- of the rris = 1/2 plateau. The Hamiltonian (|80|) 
describes an isotropic spin-1/2 antiferromagnet. From the comments made earlier, we see 
that this model only has the two saturation plateaus at rris = =tl/2, and no other plateau 
in between. For h = 0, the two-spin correlations decay as power laws with the exponent 
?7 = 1; see Eq. (pSf). 

We now use the LEH approach to study a two-chain spin-1/2 ladder with the following 
Hamiltonian, 



H = J Z_^ Sl,n ■ S2,n + "^2 2^ 2^ Sa,n " Sa,n+1 

n a=l n 

2 

+ 2-^1 zl 'S'l,n ■ 5'2,„+l - h 2_, zl ^a,n ) (§3) 



a=l n 



as shown in Fig. 9. The model may be viewed as a single chain with an alternation in 
nearest-neighbor couplings J' and 2Ji (dimerization), and a next-nearest-neighbor cou- 
pling J2 (frustration) . Eq. (^) has been studied from the point of view of magnetization 
plateaus using a first-order LEH, bosonization and exact diagonalization [|TT|, |r^, Q. 

We begin by setting Ji = J2 = 0, and studying the four states on each rung. These 
are specified by giving the configurations ± of the spins on chains 1 and 2 as follows. The 
three triplet states with S = 1 are denoted as |1), |2) and |3), where |1) = | -f -f-) and the 
other two states are obtained by acting on it successively with S~ . These three states have 
energy J'/4 in the absence of a magnetic field. The singlet state |4) = [\-\ — ) — | — \-)]/^/2 
has energy — 3J'/4. The states |1) and |4) become degenerate at a field h^ = J'. We now 
develop perturbation theory by assuming that Ji, J2 and h — h^ are all much less than J' . 
The perturbation is V^ = Y,n Ki,n+i where 

2 ^ ^ 

Vn,n+1 = J2 2^ «." ' "^a.n+l '^ ^Jl Oi^n ' 'J2,n+1 

a=l 
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Figure 9: Schematic picture of the two-chain ladder described in Eq. (^3|). The labels 1 
and 2 denote the two chains. 






M) 



The actions of this operator on the four low-energy states of a pair of neighboring rungs 
can be easily obtained. We now introduce effective spin-1/2 operators Sn on each rung 
which act on the two low-energy states. The LEH is then found to be 



H, 



eff 



J2- Jl) zK'^n'^n+l + 'S'^'S'^+l) + -j{ J2 + Jl) zl ^nSn+1 



- ih - J' 



Jl J2 ^ ST Qz 



^5) 



We now compute the field hi above which the state |111 ■ ■ ■) becomes the ground state. 
The dispersion of a spin wave, in which one rung is equal to |4) and all the others are 
equal to |1), is given by 



u{k) = h — J' — Jl — J2 + {J2 — Jl) cosk . 



m) 



By minimizing this as a function of k in various regions in the parameter space (Ji, J2), 
and then setting that minimum value equal to zero, we find that hi is given by 



hi = J' + 2Ji if J2 < Jl , 
= J' + 2J2 if J2>Ji- 



^7) 



This is the lower critical field he- of the saturation plateau with magnetization rris = 1 
per rung. Similarly, we can find the field /12 from the dispersion of a spin wave in which 
one rung is equal to |1) and the rest are equal to |4). The dispersion is given by 

uj{k) = - h + J' + {J2-J1) cosk. (88) 



By setting the minimum of this equal to zero, we find that /i2 is given by 

h2 = J' + J2 - Jl if J2 < Jl , 
= J' — J2 + Jl if J2 > Jl • 



^9) 
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This is the upper critical field hc+ of the saturation plateau with magnetization rng = 
per rung. 

Finally, we can see that the first-order terms in ( P5| ) are of the same form as the XX Z 
model in (0). We can always make the coefficient of the first term in (|85|) positive, if 
necessary by performing a rotation S^ -^ (— 1)"S'^, S^ -^ {—1)^S^ and S^ -^ S^. We 
then get a Hamiltonian of the form 

n n 

-(h - J' - J - j)T.S;.- (90) 

^ = ,f/^'n ■ (91) 

Z I J2 — "-'1 1 

From the earlier comments, we see that the two-chain ladder has an additional plateau 
at TTis = 1/2 for A > 1, i.e., if J2 + Ji > 2| J2 — Ji|. In particular, A = 00 for J2 = Ji; 
the rUg = 1/2 plateau should then extend all the way from the upper critical field of the 
rUs = plateau to the lower critical field of the m^ = 1 plateau. This can be seen in Fig. 
10 which is taken from Ref. |T^; the dimerization parameter a in that figure is related 



This is an XX Z model with 



to our couplings by J' = 1 + a and 2Ji = 1 — a. Note that the m^ = 1/2 plateau is 
particularly broad at a = 0.6, i.e., J2 = Ji = 0.2, and that it actually touches the rris = 1 
plateau on the right. The fact that it does not extend all the way up to the rris = 
plateau on the left is probably because we have ignored higher-order terms which lead to 
deviations from the XX Z model. 

To summarize, we studied a three-chain spin-1/2 ladder with a large ratio of interchain 
coupling to intrachain coupling using a LEH approach. We found a wide plateau with 
rung magnetization given by rus = 1/2. The two-spin correlations are extremely short- 
ranged in the plateau. All these are consistent with the large magnetic gap. At other 
values of m, the two-spin correlations fall off as power laws; the exponents can be found 
by using the first-order LEH which takes the form of an XX Z model in a longitudinal 
magnetic field. We also used the LEH approach to study a two-chain ladder with an 
additional diagonal interaction. In addition to a plateau at rUs = 0, this system also has 
a plateau at mg = 1/2 for certain regions in parameter space. The rris = 1/2 plateau 
is interesting because it corresponds to degenerate ground states which spontaneously 
break the translation invariance of the Hamiltonian. This can be understood from the 
LEH which, at first-order, is an XX Z model with A > 1. 



6 Summary 

We have presented some field theoretic methods for studying the properties of quantum 
spin systems in one dimension. Each of these methods has a particular regime of validity 
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Figure 10: Magnetization plateaus of the two-chain ladder as a function h and a for 
J2 = 0.2. The numbers 0, 1/2 and 1 correspond to the values of vris at the plateaus. 

(i.e., large S for the NLSMs, small S for bosonization, and weak perturbations for the 
LEHs) within which the method can give a reasonable qualitative picture of the ground 
state and low-energy excitations. Such a picture is very useful for gaining a quick under- 
standing of a given model, even though one may then need to use numerical methods like 
the DMRG to obtain quantitative results. 
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